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ABSTRACT 

We present galaxy-galaxy lensing measurements over scales 0.025 to 10 h^^ Mpc in the Sloan Digital 
Sky Survey. Using a flux-limited sample of 127,001 lens galaxies with spectroscopic redshifts and mean 
luminosity (L) ~ and 9,020,388 source galaxies with photometric redshifts, we invert the lensing 
signal to obtain the galaxy-mass correlation function ■ We find is consistent with a power-law, 
= (r/ro)-^, with best-fit parameters 7 = 1.79 ± 0.06 and tq = (5.4 ± 0.7)(0.27/17,„)i/'>'/i-i Mpc. 
At fixed separation, the ratio £,gg/S,gm — b/r where b is the bias and r is the correlation coefficient. 
Comparing to the galaxy auto-correlation function for a similarly selected sample of SDSS galaxies, 
we find that b/r is approximately scale independent over scales 0.2 — 6.7h~^ Mpc, with mean (b/r) = 
(1.3±0.2)(rim/0.27). We also find no scale dependence in b/r for a volume limited sample of luminous 
galaxies (-23.0 < Mr < -21.5). The mean b/r for this sample is {b/r)vum = (2.0 ± 0.7)(f^m/0.27). 
We split the lens galaxy sample into subsets based on luminosity, color, spectral type, and velocity 
dispersion, and see clear trends of the lensing signal with each of these parameters. The amplitude 
and logarithmic slope of S^gm increases with galaxy luminosity. For high luminosities (L ^ 5L*), ^gm 
deviates significantly from a power law. These trends with luminosity also appear in the subsample 
of red galaxies, which are more strongly clustered than blue galaxies. 

Subject headings: cosmology: observations — dark matter — gravitational lensing — large-scale struc- 
ture of the universe 



1. INTRODUCTION 

The measurement of galaxy clustering has long been a 
primary tool in constraining structure formation models 
and cosmology. Yet the power of galaxy surveys to dis- 
criminate between models is partially compromised by 
the fact that they provide only an indirect measure of 
the underlying mass distribution, subject to consider- 
able uncertainties in the bias, that is, in how luminous 
galaxies trace the mass. In this context, the galaxy-mass 
cross-correlation function can provide important ad- 
ditional information, since it is in some sense a 'step 
closer' to the clustering of mass. Moreover, comparing 
igm with the galaxy auto-correlation function yields 
a measure of the bias and therefore a constraint on the- 
ories of galaxy formation. In this paper we measure the 
correlation between galaxies and mass using weak gravi- 
tational lensing. 
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In the current paradigm of structure formation, the for- 
mation of galaxies is heuristically divided into two parts. 
On large scales, cosmological parameters and the prop- 
erties of the dark matter determine the growth of den- 
sity perturbations and the eventual formation of massive 
dark halos. On smaller scales, hydrodynamic and other 
processes shape how luminous galaxies form within dark 
matter halos and how they evolve as halos accrete and 
merge. A natural consequence of this picture is that 
the galaxy distribution is related to but differs in de- 
tail from the mass distribution. This difference arises in 
part because halos are more strongly clustered than the 
dark matter as a whole, and more massive halos are more 
strongly clustered than less massive ones (Kaiser 1984). 
In addition, the efficiency of forming luminous galaxies 
of different types and luminosities varies (primarily) with 
halo mass. 

The galaxy bias parameter can be defined in a number 
of ways, but a traditional one is as the ratio of the galaxy 
and mass auto-correlation functions at fixed separation 
(Kaiser 1984), 



gg 



(1) 



The amplitude of the galaxy-mass cross-correlation ^gm 
relative to ^gg and can be expressed in terms of the 
correlation coefficient (Pen 1998), 



gm 



(2) 



so that ^g,„ = br^mm- In general, b and r can be time- 
dependent functions of the pair separation and depend 
on galaxy properties. Since we will be comparing the 
galaxy-galaxy and galaxy-mass correlation functions, we 
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will constrain the ratio 

Galaxy surveys have provided a wealth of informa- 
tion on the behavior of the galaxy bias 6 as a function 
of galaxy luminosity and type. For example, on scales 
r < 20h~^ Mpc, the clustering amplitude ^gg increases 
with luminosity (Norberg et al. 2001; Zehavi et al. 2002; 
Norberg et al. 2002; Zehavi et al. 2004), while the am- 
plitude and shape of S^gg vary systematically from early 
to late galaxy types (Davis & Gellcr 1976; Norberg et al. 
2001; Zehavi et al. 2002). Moreover, the nearly power-law 
behavior of ^gg as well as the small departures therefrom 
(Zehavi et al. 2003), combined with the assumption that 
the dark matter distribution is described by cold dark 
matter models, indicate that the bias is scale-dependent 
on these scales. 

On larger scales, r > 20h~^ Mpc, there is evi- 
dence from higher-order galaxy correlations (Ftieman & 
Gaztanaga 1999; Szapudi et al. 2002; Verde et al. 2002), 
from the cosmic shear weak Icnsing power spectrum 
(Hoekstra et al. 2002b; Jarvis et al. 2003), and from com- 
parison of the 2dF galaxy power spectrum (Percival ct al. 
2001) with the WMAP cosmic microwave background 
(CMB) temperature angular power spectrum (Spcrgel 
ct al. 2003) that the linear galaxy bias parameter is 
of order unity for optically selected galaxies. (Here, 
b'l^^ is the ratio of the galaxy correlation function to 
the mass correlation function computed in linear per- 
turbation theory; it is related but not identical to the 
bias defined in Equation 1.) Measurement of the pa- 
rameter /3 = ^^j^ /hun from redshift space distortions 
in galaxy surveys, combined with independent evidence 
that ~ 0.3, also indicates hnn{L*) ~ 1 on large scales 
(Peacock et al. 2001). 

In recent years, weak gravitational lensing has become 
a powerful tool for probing the distribution and cluster- 
ing of mass in the Universe. We focus on galaxy-galaxy 
lensing, the distortion induced in the images of back- 
ground (source) galaxies by foreground lens galaxies. Al- 
though the typical distortion induced by a galaxy lens is 
tiny (~ 10"'^) compared to the intrinsic ellipticities of the 
source galaxies (~ 0.3), the signals from a large sample of 
lens galaxies can be stacked, providing a mean measure- 
ment with high signal to noise. The mean lensing signal 
can be used to infer the galaxy-mass cross-correlation 
function which, when compared with the galaxy auto- 
correlation function, constrains the amplitude and scale 
dependence of the bias. 

The first detection of galaxy-galaxy lensing was made 
by Brainerd et al. (1996), and the field has progressed 
rapidly since then (dell' Antonio & Tyson 1996; Griffiths 
et al. 1996; Hudson et al. 1998; Fischer et al. 2000; Wilson 
et al. 2001; Smith et al. 2001; McKay et al. 2002; Hoek- 
stra et al. 2003a). The first high S/N measurements were 
made in the Sloan Digital Sky Survey (SDSS) (Fischer 
et al. 2000) . Recent studies have benefited from improved 
data analysis and reduction techniques and from surveys 
which are specifically designed for lensing (Hoekstra et al. 
2001). Most work in galaxy lensing has concentrated on 
its power to constrain galaxy halo parameters (Brain- 
erd et al. 1996; Hudson et al. 1998; Fischer et al. 2000; 
Hoekstra et al. 2003a). However, galaxy lensing has also 



been used to measure the bias directly (Hoekstra et al. 
2001, 2002a); their results indicate that h and r are scale- 
dependent over scales ~ 0.1 — bh~^ Mpc, but that the 
ratio b/r is nearly constant at h/r ~ 1.1 over this range. 

A significant step forward in galaxy-galaxy lensing 
came with the use of samples of lens galaxies with spec- 
troscopic redshifts (Smith et al. 2001; McKay et al. 2002). 
Lensing measurements could then be made as a function 
of physical rather than angular separation, placing lens- 
ing correlation measurements on a par with the auto- 
correlation measurements from galaxy redshift surveys. 
Incorporation of photometric redshifts for the source 
galaxies (Hudson et al. 1998) also substantially reduces 
errors in the lens mass calibration due to the breadth of 
the source galaxy redshift distribution. 

In this paper, we study galaxy-mass correlations in the 
SDSS using weak gravitational lensing. Using a sam- 
ple of 127,001 galaxies with spectroscopic redshifts and 
9,020,388 galaxies with photometric redshifts, we mea- 
sure the lensing signal with high S/N over scales from 
0.025 - 10/z~^ Mpc. This is the first galaxy lensing 
study to incorporate both spectroscopic lens redshifts 
and photometric source redshifts, it is by far the largest 
galaxy lens-source sample compiled, and it extends to 
scales larger than previous galaxy-galaxy measurements. 
A similar spectroscopic sample has been used for galaxy 
auto-correlation measurements in the SDSS (Zehavi et al. 
2003). We compare the galaxy- mass and galaxy-galaxy 
correlations to constrain b/r over scales from 200 kpc to 
10 Mpc. We also use the spectroscopic and photomet- 
ric data from the SDSS to divide the lens galaxy sample 
by luminosity, color, spectral type, and velocity disper- 
sion. We see clear dependences of ^g„i on each of these 
properties. 

The layout of the paper is as follows: In §2 we intro- 
duce lensing and the measurement methods. In §3 we 

discuss the SDSS data, reductions, and sample selection. 
The basic measurement of the lensing signal AS is pre- 
sentcid in §4, and important checks and corrections for 
systematic errors using random points are discussed in 
§4.1. In §5-7 we use the data to infer the galaxy- mass 
correlation function S,gm and compare it with indepen- 
dent measurements of S,gg to constrain the bias. In §8-9 
we explore the dependence of galaxy-mass correlations 
on galaxy luminosity and type. We conclude in §10 and 
discuss possible systematic errors in the Appendix. 

Throughout this paper, where necessary we use a 
Friedman-Robertson- Walker cosmology with U,m = 0.27, 
JIa — 0.73, and Hq = 100 h km/s. All distances, densi- 
ties, and luminosities are expressed in comoving coordi- 
nates. 

2. LENSING AND GALAXY-MASS CORRELATIONS 

2.1. Gravitational Shear and the Galaxy-mass 

Correlation Function 

In this section we review the relation between the in- 
duced shear, which can be estimated from source galaxy 
shape measurements, the galaxy-mass cross correlation 
function £^gm, and the projected cross-correlation func- 
tion Wgm- The tangential shear, 7^, azimuthally aver- 
aged over a thin annulus at projected radius R from a 
lens galaxy, is directly related to the projected surface 
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mass density of the lens within the aperture, 

7T X T,„u = E(< R) - S(ii) = AS , (4) 
where E(< R) is the moan surface density within radius 
R, and is the azimuthally averaged surface density 
at radius R (Miralda-Escude 1991; Kaiser et al. 1994; 
Wilson et al. 2001). The proportionality constant T,crit 
encodes the geometry of the lens-source system, 

^crit — c^Ds ' 

where D^, Ds, and D^s are angular diameter distances 
to lens, source, and between lens and source. 

Due to the subtraction in equation 4, uniform mass 
sheets (such as the mean density of the universe p = 
^ImPcrii) do not Contribute to AS — it measures the mean 
excess projected mass density. The mean excess mass 
density at radius r from a galaxy is p £,gm(r). The mean 
excess projected density S(i?) is given by the radial in- 
tegral: 



(S(i?)) = / p(gmix,y,z)dz = pWgm{R) 



(6) 



where 



Ugm is the projected galaxy-mass correlation func- 
tion and R — {x^ + 2/^)^^^ is the projected radius. The 
observable AS is itself an integral over S(ii) and hence 

(■R 

R'dR' Wgm{R')-Wgm{R) 



(AS(i?)) 



2 







(7) 

If the cross-correlation function can be approximated 

by a power- law in separation, = (r/ro )"'•', then Wgm 
can be written as 



Wgm{R) = F{j,ro)R^- 



(8) 

where F(7,ro) = r;j'r(0.5)r[0.5(7 - l)]/r(0.57) (Davis 
& Peebles 1983). In that case, the mean lensing signal 
AS is also a power law with index 7 — 1 and is simply 
proportional to pWgm, 

{AJ:{R))=(^l-^^pWgm{R). (9) 

More generally, the three-dimensional galaxy-mass cor- 
relation function can be obtained by inverting AS di- 
rectly. Differentiating equation 7, we find 



dw. 



-p- 

The derivative dwo 



gm 



dR 



dAY, AS 
dR "^^~R~' 



(10) 



ugm/dR can be integrated to obtain ^gm 
using an Abell formula (Saunders et al. 1992): 

^.r.ir)='-rdR^,P4^ (11) 



'(i?2 -r2)i/2 

In practice, the data only cover a finite range of scales 
up to r = Rmax- The estimated £,gm, integrating to Rmax 
is related to the true ^gm by 



f''^*(r) 



gm 



TT 



dR 



-dwgnJdR 



(12) 



where the last term reminds us of the (unknown) con- 
tribution from scales beyond those for which we have 
measurements. Provided the integrand falls sufficiently 
fast with separation, this term is negligible for scales r 
smaller than a fraction of Rmax- Furthermore, since ^gm 
is linear in AS, the covariance matrix of the latter can 
be straightforwardly propagated to that of the former. 



2.2. Estimating AS 

We estimate the shear by measuring the tangential 

component of the source galaxy ellipticity relative to the 
the lens center, e+, also known as the E-mode. In gen- 
eral, the shear is related in a complex way to e+ (Schnei- 
der & Seitz 1995), but in the weak lensing regime the 
relationship is linear: 



(13) 



where e™* is the intrinsic ellipticity of the source, jt is 
the shear, and TZ is the "responsivity" (see equation 20). 
The assumption behind weak lensing measurements is 
that the source galaxies are randomly oriented in the ab- 
sence of lensing, in which case their intrinsic shapes con- 
stitute a large but random source of error on the shear 
measurement. This "shape noise" is the dominant source 
of noise for most weak lensing measurements. We discuss 
limits on intrinsic correlations between galaxy elliptici- 
ties in the appendix. 

The other component of the ellipticity, Cx, also known 
as the B-mode, is measured at 45° with respect to the 
tangent. The average B-mode should be zero if the in- 
duced shear is due only to gravitational lensing (Kaiser 
1995; Luppino & Kaiser 1997). This provides an impor- 
tant test for systematic errors, such as uncorrected PSF 
smearing, since they generally contribute to both the E- 
and B-modes. 

In order to estimate AS from the shear, we must know 
the angular diameter distances Dl, Ds, D^s for each 
lens-source pair (sec equations 4 and 5). In the SDSS, 
we have spectroscopic redshifts for all the lens galaxies, 
so that -Dl is measured to high precision (assuming a 
cosmological model). For the source galaxies, we have 
photometric redshift estimates (photoz), with typical rel- 
ative errors of 20-30% (see §3.1.5), so there is comparable 
uncertainty in the value of Y^crit for each lens-source pair. 

Given the known redshifts of the lenses, the distribu- 
tion of errors in the source galaxy ellipticity, and the dis- 
tribution of errors in the photometric redshift for each 
source, we can write the likelihood for AS from all lens- 
source pairs. 



£(AS) 



n 



n / 



dzip{zi)p{^i.\ziz{) , (14) 



where 7^ = e\/2TZ is the shear estimator for the «th 
source galaxy, P{zl) is the probability distribution for its 
redshift (the product of the Gaussian error distribution 
returned by the photoz estimator and a prior based on 
the redshift distribution for the source population; see 

§3.1.5), and f (iTkli-^i) is the probability distribution 
of the shear given the source and lens redshifts, which is 
a function of the desired quantity AS: 



4) ^ 




-ASxS,-4(4,4) 



(15) 



In equation 15, 4 x 0-2(7^) = a'^{e\)+a'g^: the shear un- 
certainty is the sum of the measurement variance a'^{e\^) 
and the intrinsic variance in the shapes of the source 
galaxies a'gj^ = ((e™*)^). The shape noise measured 
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from bright, well-resolvcd galaxies is ctsn ~ 0.32, and 
the typical measurement error cr(e+) ranges from ~0.05 
for r=18 to ~0.4 for r=21.5. The intrinsic shape dis- 
tribution is not Gaussian as we have assumed in equa- 
tion 15, but it is symmetric. Monte Carlo simulations 
indicate that this approximation does not bias the mea- 
surement of AS within our measurement uncertainties, 
provided we take asN as the standard deviation of the 
non-Gaussian shape distribution. 

Although the typical uncertainty in the source galaxy 
photometric redshifts is 20-30%, this is small compared 
to the relative shear noise, which is typically cfsnHt ^ 
300%. Assuming shape error is the dominant source of 
noise, we can approximate equation 14 as 



log£(AS) = ^ - 



(16) 



where we now use the critical surface density averaged 
over the photoz distribution for each source galaxy, 

{Kru)j,i = j d4P{zi)^-J^{zlz{) (17) 

Monte Carlo simulations indicate that the true likelihood 
approaches this Gaussian approximation after stacking 
only a few hundred lenses. 

The maximum likelihood solution is the standard 
weighted average, 



AE 



W 



(18) 



where 



AS,.,=75^/(S,-4),., 



(19) 



Although this simple inverse variance weighting is not 
optimal (Bernstein & Jarvis 2002), it does lead to unbi- 
ased results and only a slight increase in the variance of 
AS. 

As indicated by equation 13, the ellipticity induced by 
a shear depends on the object's shear responsivity TZ. 
A measure of how an applied shear alters the shape of 
a source, TZ depends on the object's intrinsic ellipticity 
and is similar to the shear polarizability of Kaiser et al. 
(1995). Following Bernstein & Jarvis (2002), we calcu- 
late a mean responsivity as a weighted average over all 
tangential ellipticities. 



T.j,iWj,i [1 -fco -^1(6;)^^] 



(20) 



where the weights are the same as in equation 19. We 
have again assumed a Gaussian distribution of elliptici- 
ties so that fco and fci are simple. 



^o = (l - /)(7siv> ki=f, 



f-- 



(21) 



The mean TZ for our sources is 0.86. 



3. DATA 

The Sloan Digital Sky Survey (SDSS; York et al. 

(2000)) is an ongoing project to map nearly 1/4 of the 
sky in the northern Galactic cap (centered at 12*^20™, 
-1-32.8°). Using a dedicated 2.5 meter telescope lo- 
cated at Apache Point Observatory in New Mexico, the 
SDSS comprises a photometric survey in 5 bandpasses 
{u, g,r,i, z; Fukugita et al. (1996)) to r '--^ 22 and a 
spectroscopic survey of galaxies, luminous red galaxies, 
quasars, stars, and other selected targets. In addition, 
the survey covers 3 long, thin stripes in the southern 
Galactic hemisphere; the central southern stripe, cover- 
ing 200 square degrees, will be imaged many times, al- 
lowing time-domain studies as well as a deeper co-added 
image. 

We select our lens galaxies from the SDSS main galaxy 
spectroscopic sample. For our background sources we use 
only well-resolved galaxies drawn from the photometric 
survey with well- measured photometric redshifts. Each 
sample is described in detail below. 

3.1. Imaging Data 

Imaging data are acquired in time-delay-and-integrate 
(TDI) or drift scan mode. An object passes across the 
camera (Gunn et al. 1998) at the same rate the CCDs 
read out, which occurs continuously during the exposure. 
The object crosses each of the 5 SDSS filters in turn, 
resulting in nearly simultaneous images in each bandpass. 
In order for the object to pass directly down the CCD 
columns, the distortion across the field of view must be 
exceedingly small. This is advantageous for lensing, since 
distortions in the optics cause a bias in galaxy shapes 
(§3.1.3). The distortion in the optics of the SDSS 2.5m 
telescope is negligible (Stoughton et al. 2002) and can be 
ignored. 

The imaging data are reduced through various soft- 
ware pipelines, including the photometric (PHDTD, Lup- 
ton et al. (2001)), astrometric (Pier et al. 2003), and cal- 
ibration (Hogg et al. 2001; Smith et al. 2002) pipehnes, 
leading to calibrated lists of detected objects. The cali- 
brated object lists are subsequently fed through various 
target selection pipelines (Eisenstein et al. 2001; Strauss 
et al. 2002; Richards et al. 2002) which select objects for 
spectroscopic followup. 

The shape measurements discussed in §3.1.1 are im- 
plemented in PHOTO (v5_3), so we work directly with the 
calibrated object lists. These lists contain, among many 
other things, position (RA,DEC), several measures of the 
flux, diagnostic flags for the processing, and moments of 
the light distribution for each object and for the local 
Point Spread Function (PSF) (Stoughton et al. 2002). 
We augment the parameters measured by PHOTO with 
the probability that each object is a galaxy (§3.1.4) and 
with photometric redshifts. 

3.1.1. Shape Measurements 

Weak lensing measurements rely on the assumption 
that source galaxy shapes arc an unbiased, albeit rather 
noisy, measure of the shear induced by foreground lenses. 
Therefore high S/N shape measurements and accurate 
corrections for bias are crucial for weak lensing measure- 
ments. For both shape measurements and corrections 
we use techniques described in Bernstein & Jarvis (2002) 
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(hereafter BJ02). 

Wc determine the apparent shapes of objects from their 
flux-weighted second moments. 



Qr. 



(22) 



where Im,n is the intensity at pixel m, n and Wm,n is an 
elliptical Gaussian weight function, iteratively adapted 
to the shape and size of the object. Initial guesses for the 
size and position of the object are taken from the Pet- 
rosian radius and PHOTO centroid. Objects are removed if 
the iteration does not converge or if the iterated centroid 
wanders too far from the PHOTO centroid. The shape is 
parametrized by the polarization, or ellipticity, compo- 
nents, defined in terms of the second moments, 



-Q 



22 



62 = 



Qii + Qii 

2Ql2 

Qii + Q22 



(23) 

(24) 



The polarization is related directly to the shear via equa- 
tion 13. 

The errors in the moments are calculated from the pho- 
ton noise under the assumption that each measurement 
is sky noise limited, which is a good approximation for 
the faint sources. Because the same pixels are used for 
measuring both components of the ellipticity, there is a 
small covariancc between the components of the ellip- 
ticity. This covariance is also calculated and properly 
transformed when rotating to the tangential frame for 
shear measurements. 

Because the natural coordinates for the SDSS are sur- 
vey coordinates (A, 77) (Stoughton et al. 2002), wc ro- 
tate the cllipticities into that coordinate system for the 
lensing measurements. The full covariance matrix for 
(61,62) is used to transform the errors under rotations. 
We combine the shape measurements from the g, r, and i 
bandpasses using the covariancc matrices. This increases 
the S/N, simplifies the analysis, and reduces bandpass- 
dependent systematic effects. The u and z bands are 
much less sensitive and would contribute little to the 
analysis. 

3.1.2. PSF Reconstruction 

An anisotropic PSF, caused by instrumental and atmo- 
spheric effects, smears and alters the shapes of galaxies 
in a way which can mimic lensing. In addition, the finite 
size of the PSF — the seeing — tends to circularize the im- 
age, reducing the measured ellipticity. In order to correct 
for these effects, one should, in principle, determine the 
exact shape and size of the PSF at the position of every 
source galaxy. Imaging in drift scan mode produces im- 
ages that are long, thin stripes on the sky. Because the 
PSF varies over time (along the direction of the scan) 
as well as across the camera, it must be tracked as a 
function of position in the overall image (see §3.1.3). 

The photometric pipeline uses Karhunen-Loevc (KL) 
decomposition (Hotelling 1933: Karhunen 1947: Loeve 
1948) to model the PSF (Lupton et al. 2001). The PSF 
is modeled on a frame- by- frame basis, where frames are 
defined as 2048 x 1490 pixel chunks composing the long 
SDSS image. A set of bright, isolated stars are chosen 
from ±2 frames around the central frame, with typically 



15 — 25 stars per frame. These stellar images P(j)(M,w) 
are used to form a set of KL basis fmictions or eigen- 
images Br{u,v), in terms of which the images can be 
reconstructed by keeping the first n terms in the expan- 
sion, 

n-l 

P^i){u,v) = J2ahBr{u,v) (25) 

r=0 

where P(i) denotes the z*'' star, and u, v arc pixel posi- 
tions relative to the object center. 

The spatial dependence of the coefficients a^^^ are de- 
termined via a polynomial fit. 



Hi) 



l+m<N 

E 

l=m=Q 



■{i)y{i) 



(26) 



where x, y are the coordinates of the center of the i star 
relative to the center of the frame, N is the highest order 
in X, y included in the expansion, and are determined 
from minimizing 



n-l 



2 



X = 



(27) 



Only the stars on ±1/2 a frame surrounding the given 

frame are used to determine the spatial variation, while 
stars from ±2 frames are used to determine the KL basis 
functions. 

In PHOTO, the number of terms used from the KL basis 
is usually n = 3; the order of the spatial fit is A'' = 2 
unless there are too few stars, in which case the fit may 
be order 1 or even (rare). To determine the coefficients 
6[^, a total of n(A''-|- l)(A'-|-2)/2 constraints are needed, 
which may seem like too many for the typical number 
of available stars. There are many pixels in each star, 
however, so the number of spatial terms (A^ + 1)(A + 
2)/2 = 6 (for quadratic ffis) should be compared with 
the number of available stars. 

The PSF is reconstructed at the position of each object 
and its second moments are measured; these are used in 
the analytic PSF correction scheme described in §3.1.3. 

3.1.3. Shape Corrections 

To correct galaxy shapes for the effects of PSF dilution 
and anisotropy, we use the techniques of BJ02 with the 
modifications specified in Hirata & Seljak (2002). Rather 
than a true deconvolution, this is an approximate ana- 
lytic technique. The PSF is modeled as a transformation 
of the preseeing shape; to remove the effects of the PSF, 
the inverse transform must be calculated. This trans- 
formation is performed in shear space (or polarization 
space): 



R * [i-epsp) ® e] = {-epsp) ® eo 



(28) 



where e is the preseeing polarization, cq is the measured 
polarization, epsp is the polarization of the PSF, and 
R is the resolution parameter, which is related to the 
polarizability of Kaiser et al. (1995). The operator © is 
the shear addition operator defined in BJ02. If the shear 
addition operator were a simple addition, then equation 
28 would reduce to: 



6 = [60 - (1 - R)epsF]/R 



(29) 
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For unweighted moments {Wm.n = 1) or if the objects 
and PSF have Gaussian surface brightness profiles, the 
formulae in BJ02 are exact, and the resolution parame- 
ter in this case is just -R = 1 — {spsF I SobjY ■, where s is 
the linear size of the object or PSF. The light profiles of 
galaxies and of the PSF differ significantly from Gaus- 
sian, however. Furthermore, as noted in §3.1.1, we use a 
Gaussian radial weight function to optimize the S/N of 
object shape measurements. As a consequence, the res- 
olution parameter R must be derived in an approximate 
way that accounts for both the weight function and the 
non-Gaussianity of the light profiles. For this work we 
use a weighted fourth moment to correct for higher order 
effects as discussed in BJ02 and Fischer et al. (2000). 

3.1.4. Star-Galaxy Separation 

To separate stars and galaxies cleanly at all magni- 
tudes, we use the Bayesian method discussed in Scranton 
et al. (2002). The method makes use of the concentra- 
tion parameter, which can be calculated from parameters 
output by PHOTO. The concentration is the difference be- 
tween the object's PSF magnitude and exponential disk 
magnitude. The PSF magnitude is derived by fitting the 
local PSF shape to the object's light profile. The only 
free parameter in this fit is the overall flux. An expo- 
nential disk is also fit to the object, but in that case 
the scale length is also a free parameter. Thus, large 
objects have more flux in the exponential than the PSF 
fit and correspondingly large concentration, while stars 
have concentration around zero. 

At bright magnitudes, galaxies and stars separate 
cleanly in concentration space. At fainter magnitudes, 
photometric errors increase and the distributions overlap. 
This is demonstrated in figure 1, which shows the distri- 
bution of concentration for objects with 20 < r < 21 
and 21 < r < 22 drawn from 100 frames of a single SDSS 
imaging run (3325), with mean seeing of 1.25" (typical of 
SDSS image quality). In bad seeing conditions, a smaller 
percentage of galaxies are larger than the PSF, again 
making separation difficult. 

If we know the distribution in concentration for galax- 
ies and stars as a function of seeing and magnitude, we 
can assign each object a probability Pg that it is a galaxy. 



500 




-0.2 0.0 0.2 0.4 0.6 

concentration 

Fig. 1. — Concentration distribution for objects with 20 < r < 
21 (dark curve) and 21 < r < 22 (light curve). Stars have 
concentration near zero. At faint magnitudes, stars and galaxies 
are not as easily separated. 



given its concentration, magnitude, and the local seeing. 
For our source galaxy sample, we can then select objects 
which have high values of Pg. To map out the distri- 
bution in concentration, we use regions from the SDSS 
Southern Survey which have been imaged many times. 
Some of these regions have been imaged as many as 16 
times, with an average of about 8. By averaging the 
flux for each object from the multiple exposures to ob- 
tain higher S/N, a clean separation between the star and 
galax;y concentration distributions is achieved to fainter 
magnitudes. We have maps of the concentration distri- 
bution for 16 < r < 22 and 0.9 < seeing < 1.8, 
allowing us to accurately calibrate Pg and therefore de- 
fine a clean galax;y sample. In the few regions of very 
bad seeing, we extrapolate the concentration conserva- 
tively, erring on the side of including fewer galaxies in 
the sample. In section 3.1.6 we discuss the cuts used for 
our source sample. 

3.1.5. Photometric Redshifts 

A photometric redshift (photoz) is estimated for each 
object in the source catalog. The repaired template fit- 
ting method is used, described in detail in Csabai et al. 
(2000) and implemented in the SDSS EDR (Csabai et al. 
2003) as well as the DRl (Abazajian et al. 2003). This 
technique uses the 5-band photometry for each object as 
a crude spectrum. The algorithm compares this spec- 
trum to templates for different galaxy types at different 
redshifts. The result is an estimate of the type and red- 
shift of each galaxy. 

There is a large covariance between the inferred type 
of the galaxy and its photometric redshift. The code 
outputs a full covariance matrix for type and redshift. 
Because we do not use the type information, we use the 
error marginalized over type. We further assume that the 
resulting error is Gaussian, which is only a good approx- 
imation for high S/N measurements. This introduces a 
bias in the estimate of (S^^^^f), but Monte Carlo simula- 
tions indicate that this is a negligible effect. 

From comparisons to galaxies with known redshifts, 
the rms in the SDSS photoz estimates is found to be 
~ 0.035 for r < 18, increasing to ~ 0.1 for r < 21. About 
30% of our source galaxy sample has r between 21 and 22. 
Although the photozs are less reliable in this magnitude 
range, these objects receive relatively little weight in the 
analysis, because they have large shape errors and large 
photoz errors (and hence small (EJ^^J). 

The photoz distribution for the sources used in this 
study is shown in figure 2. The histogram shows the 
photometric redshifts, and the smooth curve is the dis- 
tribution calculated by summing the Gaussian distribu- 
tions for each galaxy. We use this smooth curve as a prior 
on the photometric redshift when calculating the inverse 
critical density for each lens-source pair (see §2.2). The 
large peaks in the distribution are most likely not real, 
but rather the result of degeneracies in the photometric 
redshift estimation. This issue is addressed further in the 
appendix. 

Although there is significant overlap between the dis- 
tribution of photozs and the distribution of lens redshifts 
(figure 4), sources with photozs in front of or near the lens 
redshift are given appropriately small weight according 
to equation 17. 
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Fig. 2. — Distribution of photometric redshifts for sources used 
in this study. The histogram shows the photozs in bins of Az = 
.01, and the smooth curve is derived from summing the Gaussian 
distributions associated with each object. 



3.1.6. Defining the Source Sample 

Source galaxies are drawn from SDSS imaging stripes 
9-15 and 27-37, covering a region of nearly 3800 square 
degrees. An Aitoff projection displaying the positions of 
these sources (as well as the lenses) is shown in figure 3. 
We make a series of cuts aimed at ensuring that the sam- 
ple is of high purity (free from stellar contamination) and 
includes only well-resolved objects with usable shape in- 
formation. We first require that the extinction-corrected 
r-band Petrosian magnitude is less than 22. We next 
make an object size cut, requiring that the resolution 
parameter R > 0.2. This removes most of the stars and 
unresolved galaxies from the sample. However, at faint 
magnitudes (r> 21), many stars and galaxies have sim- 
ilar values of R due to measurement error, so a further 
cut is needed. We employ the Bayesian galaxy proba- 
bility (§3.1.4) and find that the combination R > 0.2, 
Pg > 0.8 guarantees that the source galaxy catalog is 
greater than 99% pure for r < 21.5 and greater than 
98% pure for 21.5 <r< 22. 



those with photoz errors greater than 0.4, and we further 
exclude objects with photoz less than 0.02 or greater than 
0.8 since failed measurements tend to pile up at a photoz 
of 0.0 or 1.0. This removes another 10% of the objects. 
The final source catalog contains 9,020,388 galaxies, cor- 
responding to a density of about 1 — 2 source galaxies 
per square arcminute, depending on the local seeing. 

3.2. Spectroscopic Data 

The lens galaxies are selected from the SDSS "main" 
galaxy spectroscopic sample, which is magnitude- 
(14.5 <r< 17.77) and surface brightness-limited {fir < 
23.5), although these limits varied during the commis- 
sioning phase of the survey. See Strauss et al. (2002) for 
a description of "main" galaxy target selection. 

SDSS spectroscopy is carried out using 640 optical 
fibers positioned in pre-drilled holes in a large metal 
plate in the focal plane. Targeted imaging regions are 
assigned spectroscopic plates by an adaptive tiling algo- 
rithm (Blanton et al. 2003), which also assigns each ob- 
ject a fiber. The spectroscopic data are reduced to 1-d 
spectra by SPECTR02d, and the SPECTROld pipchne out- 
puts redshift and an associated confidence level, spectral 
classification (galaxy, quasar, star), line measurements, 
and spectral type for galaxies, among other parameters 
(Stoughton et al. 2002). In addition, the velocity dis- 
persion is measured for a large fraction of the early type 
galaxies. 

For this analysis, we use a subset of the available spec- 
troscopic "main" galaxy sample known as LSS samplel2 
(M. Blanton 2003, private communication). Although 
we draw from a larger sample, the mask (see §3.2.3) 
was produced for this subset. This sample is also be- 
ing used for analysis of the galaxy auto-correlation func- 
tion (Zehavi et al. 2004), while a slightly earlier sample 
(scLmplell) has been used to estimate the galaxy power 
spectrum (Tegmark et al. 2004) . Using this sample allows 
us to make meaningful comparisons between the auto- 
correlation function and the galaxy-mass cross correla- 
tion function. The spectroscopic reductions used here are 
those of the SPECTROld pipeline, for which redshifts and 
spectroscopic classifications differ negligibly from those 
in the above references. 




Fig. 3. — Aitoff projection showing the positions of the sources 
(black) and lenses (red) used in this study. The section covering 
the equator is stripes 9-15. The higher latitude section is stripes 
27-37. 

Additionally, we remove about ~ 8% of the sources — 



3.2.1. Redshifts 

The SDSS spectra cover the wavelength range 3800- 
9200A with a resolving power of 1800 (Stoughton et al. 
2002). Repeated 15-minute exposures (totaling at least 
45 minutes) are taken until the cumulative median 
(S/N)^ per pixel in a fiber aperture is greater than 15 
at 5 = 20.2 and i = 19.9 in all 4 spectrograph cameras. 
Redshifts are extracted with a success rate greater than 
99%, and redshift confidence levels are greater than 98% 
for 95% of the galaxies. Repeat exposures of a number 
of spectroscopic plates indicates that "main" galaxy red- 
shifts are reproducible to 30 km/s. We apply a cut on 
the redshift confidence level at > 75%, which removes 
0.5% of the galaxies. Many galaxies are further removed 
during the lensing analysis, as discussed in §3.2.3 and 
§3.2.4. A redshift histogram is shown in figure 4 for the 
remaining 127,001 galaxies used in this work. The mean 
redshift for our sample using the relative weights from 
the lensing analysis is (z) — 0.1. 



8 



Sheldon et al. 




Fig. 4. — Redshift distribution for galaxies used in this study. 
This sample contains only "main" galaxy targets. 



3.2.2. K-corrections 

We apply K-corrections using the method discussed in 
Blanton et al. (2002) (kcorrect vl_10). Linear combi- 
nations of four spectral templates are fit to the five SDSS 
magnitudes for each galaxy given its redshift. Rest-frame 
absolute magnitudes and colors are then calculated. Fig- 
ure 5 shows the distribution of absolute Petrosian mag- 
nitude for each of the 5 SDSS bandpasses. 

3.2.3. Spectroscopic Masks 

Because SDSS spectroscopy is taken through circular 
plates with a finite number of fibers of finite angular size, 
the spectroscopic completeness varies across the survey 
area. The resulting spectroscopic mask is represented by 
a combination of disks and spherical polygons (Tegmark 
et al. 2004). Our spherical polygon mask contains 3844 
polygons covering an area of 2818 square degrees. Each 
polygon also contains the completeness, a number be- 
tween and 1 based on the fraction of targeted galaxies 
in that region which were observed. We apply this mask 
to the spectroscopy and include only galaxies from re- 
gions where the completeness is greater than 90%. The 
same criteria are used to generate the random points as 
discussed in 4.1. 




-18 -20 -22 
M - 5 log(h) 

Fig. 5. — Distribution of rest-frame absolute magnitude in each 
of the 5 SDSS bandpasses for lenses used in this study. The filled 
circles mark the magnitude range for the flux limited sample of 
Zehavi et al. (2003). Ninety percent of our lens sample falls within 
this range. 



3.2.4. Photometric Masks 

Although the correction for PSF smearing (§3.1.3) sub- 
stantially reduces the galaxy shape bias, it does not 
completely eliminate it — generally, a small, slowly vary- 
ing residual remains. Fortunately, a residual PSF bias 
that is constant over the lens aperture cancels on aver- 
age from the azimuthally averaged tangential shear, since 
two aligned sources separated by 90° relative to the lens 
contribute equally but with opposite sign. To take ad- 
vantage of this cancellation, we divide the source galaxies 
around each lens into quadrants and demand that at least 
two adjacent quadrants are free of edges and holes out 
to the maximum search radius. 

We represent the geometry of the sources using the hi- 
erarchical pixel scheme SDSSPix^ (Scranton et al. 2003), 
modeled after a similar scheme developed for CMB anal- 
ysis (Gorski et al. 1998). This scheme represents well the 
rectangular geometry of the SDSS stripes. The mask is a 
collection of pixels at varying resolution covering regions 
with holes or edges. We do not mask out bright stars 
(which cover only a tiny fraction of the survey area) for 
the lensing analysis. 

After checking the spherical polygon masks (§3.2.3), 
each lens galaxy is checked against the pixel mask to 
guarantee that it is within the allowed region. Each 
quadrant around the lens is then checked to determine 
if it contains a hole or an edge. Lenses are excluded 
from the sample if there are no adjacent quadrants that 
are completely unmasked. In addition to this cut, we 
demand that the angular distribution of source galaxies 
around the lens have ellipticity no greater than 20% in 
order to ensure the availability of pairs with 90° separa- 
tion. The same criteria are applied to the random points 
(see §4.1). 

We draw the final spectroscopic data set from SDSS 
stripes 9 — 12 (near the equator) and 28 — 37. We do 
not use galaxies from the 3 southern stripes (76,82,86), 
because they contain few lens-source pairs at large sep- 
aration. Figure 3 shows the distribution of the lens and 
source galaxies as an Aitoff projection. After applying 
the cuts described above, the final lens sample contains 
127,001 galaxies. 

4. RESULTS: THE MEAN AS FROM .02 - lOh-^ Mpc 

The mean lensing signal AS for the full sample of SDSS 
lens galaxies is shown in figure 6. The signal is unam- 
biguously detected from 25 kpc to 10 /i^^Mpc. The 
corresponding mean shear "fx is shown on the right axis. 
Corrections have been made to this profile as described 
in §4.1. These corrections are relatively small in all radial 
bins. 

The errors in figure 6 come from jackknife re-sampling. 
Although we expect statistical errors to dominate over 
sample variance even on the largest scales shown here, 
there are in addition large variations in the systematic 
errors. Due to gaps between the 5 CCD columns, two 
interleaving imaging runs make up a contiguous imag- 
ing stripe, and they are generally taken on different 
nights under different photometric conditions. As a re- 
sult, the residuals from the PSF correction vary between 
the columns of interleaving runs, each of which is ~ 0.2 

^ http:/ /lahmu. phyast.pitt.edu/~scranton/SDSSPix 
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Fig. 6.— Mean AS = E{< R) - measured for the full lens 

sample. The solid line is the best-fitting power law ASoc R~^-^^. 
The right axis shows the corresponding tangential shear 77'- 



degrees wide. The residuals also vary over time along 
the direction of the scan, with a typical scale of a few 
degrees. Thus the proper subsample size to account 
for this variation is about a square degree. We divide 
the sample into 2,000 disjoint subsamplcs, each approx- 
imately a square degree in size (see §3.2.3) and remea- 
sure AE 2,000 times, leaving out each subsample in turn. 
Figure 7 shows an image of the resulting dimension- 
less correlation matrix Corrij = Vi j / ^JVijyj~j , where 
K,, = ((AE(i?0 - (AS(i?.)'))(AS('i?,) -'(AS(i?,)))). 
The off-diagonal terms are negligible in the inner bins 
but become important beyond R ~ lh~^ Mpc. The full 
covariance matrix is used for all model fitting. 

The mean AS for the full sample is well described by 
a power law 

AS(J?) = A(i?/lMpc)-« (30) 

with a = 0.76 ± 0.05 and A = (3.8 ± O.4)/iM0pc-2. The 
outliers at intermediate radii make the reduced for 
the power law fit somewhat poor but not unacceptable: 




-1.5 -1.0 -0.5 0.0 0.5 
log(R [r' Mpc]) 

Fig. 7. — Correlation matrix for AS in figure 6 calculated using 
jackknife re-sampling. 



X^/i^ = 1.26. There is a 20% chance of exceeding this 
value randomly. 

An important check for systematic errors is the "B- 
mode", ASx the average shear signal measured at 45° 
with respect to the tangential component. If the tan- 
gential shear signal is due solely to lensing, the B-mode 
should be zero, whereas systematic errors generally con- 
tribute to both the E- and B-modes. The top panel of 
figure 8 shows AEx measured using the same source- lens 
sample used in figure 6. This measurement is consistent 
with zero. 

4.1. Systematics Tests with Random Points 

By replacing the lens galaxies with sets of random 
points, we can gauge two systematic effects on AS: resid- 
uals in the PSF correction and the radial bias due to 
clustering of sources with the lenses. We generated a 
random sample with ten times as many points as the 
lens sample, using the same masks and selection crite- 
ria described in §3.2.3 and §3.2.4. The random points 
are drawn from the same redshift distribution as the lens 
galaxies. These criteria guarantee that the same regions, 
and thus roughly the same systematics, are sampled by 
the lenses and the random points. Any non-zero lensing 
signal for the random points we ascribe to residuals in 
the PSF correction. 

The bottom panel of figure 8 shows AS measured 
around 1,270,010 random points, with errors from jack- 
knife re-sampling. Note that the random sample is large 
enough that in this case sample variance dominates the 
error in the outer bins. There is a significant signal at 
large radius. The signal at smaller radii is less well de- 
termined, but it is in any case far below the signal due 
to lenses. Interpreting the large-scale signal as system- 
atic error, we subtract it from the AS measured around 
lenses and add the errors aj^^^l = al^^ + a^^^^; this 
correction is incorporated in figure 6. 
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Fig. 8. Two tests for systematics in the lensing measurement. 
The top panel shows the "B-niode" for lensing, ASx, measured 
around the same lenses used in figure 6. The measurement is con- 
sistent with zero, as expected for lensing. The bottom panel shows 
AS measured around 1,270,010 random points. The detection at 
large radius is indicative of systematic errors, most likely from 
residuals in the PSF correction. This has been subtracted from 
the signal around lenses for figure 6. 



10 



Sheldon et al. 



The second systematic probed by the random sample 
involves clustering of the source galaxies with the lenses. 
The calculation of the mean inverse critical density in 
§2.2 properly corrects for the fact that a fraction of the 
source galaxies are in front of the lenses, but only un- 
der the assumption that the lens and source galaxies are 
homogeneously distributed. Since galaxies are clustered, 
a small fraction of the sources are in fact physically as- 
sociated with the lenses, causing a scale-dependent bias 
of the lensing signal. We correct for this by estimating 
the excess of sources around lenses compared with the 
random points. The correction factor is the ratio of the 
sums of the weights for sources around lenses and around 
random points: 



C{R) 



Nlens Y^k,l 



(31) 



where i,j indicates sources found around lenses, k,l in- 
dicates sources found around random points, and Wi_j is 
the weight for the lens-source pair (see §2.2). The value 
of C(i?) — 1 is shown for the full sample in figure 9. The 
use of photometric redshifts reduces the correction signif- 
icantly, since sources associated with the lenses are given 
little weight. We find that C{R) — 1 falls roughly as 
^-o.9_ rpj-^g correction for clustered sources is essentially 
negligible beyond 50 kpc. The signal in figure 6 has been 
multiplied by C{R). 

We measure this correction factor separately for each 
of the lens subsamples presented in later sections. Since 
galaxy clustering increases with luminosity, and higher 
luminosity lens galaxies are seen out to higher redshift 
where more of the faint sources are near the lenses, the 
correction factor C{R) increases with the luminosity of 
the lenses. The correction for our highest luminosity 
samples is 2.0 at 25 kpc compared to 1.1 for the full sam- 
ple. Similarly, the correction for early type lens galaxies 
is larger, while the correction for late types is smaller 
than for the full sample. Although the corrections for 
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Fig. 9. — Correction factor for the clustering of sources around 
the lens galaxies. The function C{R) — 1 is essentially a weighted 
cross-correlation function between lenses and sources. The AS in 
figure 6 has been multiplied by C{R), which is a negligible correc- 
tion for radii larger than ~ 50 kpc. 



some samples arc large, the value of the correction is 
well measured in each case. 

We have also calculated the average AS for the lumi- 
nosity subsamples defined below (§8). Because C{R) de- 
pends on luminosity, the average of the subsamples could 
in principle differ from the mean AE estimated from the 
full sample, as suggested by Guzik & Seljak (2002) re- 
garding the results of McKay et al. (2002). However, we 
find no significant difference in the two methods for our 
weighting scheme, which results in smaller corrections 
than those in McKay et al. (2002). 

Note our calculation of C{R) may be a slight over- or 
under-estimate, since lensing itself may induce some cor- 
relation between unassociated background sources and 
the positions of the foreground lenses through magni- 
fication. Lensing magnification can bring galaxies into 
a magnitude limited sample that would otherwise have 
been too faint to be included, and this effect will be a 
function of scale. On the other hand, the geometric dis- 
tortion induced by the lens moves the apparent position 
of background galaxies radially outward, decreasing the 
number density. The net change in number density de- 
pends on magnification /i and the slope of the galaxy 
number counts s: the ratio of counts with and with out 
lensing is N/Nq oc ^(2-5s-i) (Broadhurst et al. 1995). 
The magnification is of the same order as the shear, 
which for our lens sample is < 10^"^, and the slope s 
is typically about 0.4, which results in negligible magni- 
fication bias. 

5. THE GALAXY-MASS CORRELATION FUNCTION 

5.1. Power-law Fits to AS 

As noted in §2.1, a power-law AS is consistent, within 
the errors, with the galaxy-mass correlation function 



igm also being a power- law, = (r/ro) 



with 



slope 7 = 1 + a. Fitting for 7 and ro, we find best 
fit marginalized values of 7 = 1.76 ± 0.05 and ro = 
(5.7±0.7)(0.27/f7™)i/^/i-i Mpc. Because AS is propor- 
tional to p, we have assumed a fiducial value for the mean 
density p = pcrit^m determined by WMAP -|- ACBAR -|- 
CBI in combination with the power spectrum from 2dF- 
GRS and Lyman a data, which yields Vtm = 0.27 ± 0.02 
(Spergel et al. 2003). Marginalizing over this small un- 
certainty in does not change our error estimates. 



5.2. Inversion to ^, 



gm 



In addition to fitting a power law to AS, we performed 
a direct inversion of AS to as outlined in §2.1. 
for the main sample is shown in Figure 10. Again, since 
our measurements scale linearly with Vt^, we have as- 
sumed a fiducial value of = 0.27. 

In principle, the inversion from AS to S^gm at a given 
radius r requires data for AS to i? = 00. In practice, 
the contribution from scales beyond our data is negligi- 
ble in all radial bins but the last few. We estimate the 
contribution from large scales by extrapolating the best- 
fit power law. Because the last point would be based 
entirely on the extrapolation, we exclude it from figure 
10. The contribution of the extrapolation is small for the 
remaining points except for the last bin at 6.7 Mpc, for 
which it is a 40% effect. 

Because we do not have any a priori knowledge of 4gm 
at large scales, we must include this extrapolation in 
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Fig. 10. — Galaxy-mass cross-correlation function ^gm, obtained 
by inverting AE. The amplitude of ^gm is degenerate with Qm, 
which we have set to 0.27. The solid line is the best-fit power law 
^gm = (r/5.4/i-i Mpc)-i '^9. 



the error budget. A lower bound on ^g„i is obtained 
by assuming there is no signal beyond Rmax 

at all: 

4r = Cgt'"°'' - ^(Csm), where a{^gm) is the error on 
C(JJi^"°^ propagated from AS. The most conservative 
error estimate is then \^gm—^^gm\- This is the error plot- 
ted in figure 10 and is used for scaling the correlation 
matrix discussed below. Again, this only makes an sig- 
nificant difference in the last radial bin. This may still 
underestimate the upper bound on if in fact ^gm is 
considerably flatter at separations larger than Rmax- 

Due to the finite binning of AS (which is well described 
by a power law), linear interpolation results in a slight 
overestimate of the derivative dAT,/dR. A small 4%) 
correction for this effect has been applied in figure 10. 

Because the inversion is a linear operation on AS, it 
is straightforward to obtain the covariance matrix for 
£^gm ■ The corresponding dimensionless correlation matrix 
is shown in figure 11. Due to the differentiation and 
integration involved in the inversion, adjacent points of 
are strongly correlated; i.e., the resulting covariance 
matrix for is less diagonal than the covariance matrix 
for AS. The full covariance matrix is used for the model 
fits discussed below. 

We fit to a power law model ^g„i = (?'/?'o) For 
the inversion to ^g„j we cannot use the last radial bin, 
so this fit only applies for R < 6.7 Mpc. The results are 
summarized in table 1. We find best-fit parameters ro = 
(5.4±0.7)(0.27/f7„)i/'' and 7 = 1.79±0.06, which are in 
agreement with the results we obtained fitting directly to 
AS. These results are also in excellent agreement with 
the auto-correlation function from Zehavi et al. (2003). 
In the next section, we compare £^gm and directly to 
constrain the bias b/r. 

6. BIAS 

The ratio of the auto-correlation function ^gg to the 
galaxy-mass correlation function £^gm is a measure of the 
bias between galaxies and mass ^gg/^gm = b/r (see §1). 



It is important to use the same galaxies for both ^gm and 
^gg. We use the auto-correlation hmction from Zehavi 
et al. (2003) (hereafter Z03), measured for a flux hmited 
sample of 118,000 SDSS "main" galaxies. The selection 
function for the Z03 sample is very similar to that of 
our lens sample, but the former covers a more restricted 
range of absolute magnitudes: —22.2 < Mr < —18.9, 
while our sample spans the range —24 < Mr < —17. 
However, figure 5 shows that 90% of the galaxies in our 
lens sample in fact lie within the Z03 magnitude range. 
Also, while each galaxy is weighted differently in the two 
correlation function measurements, the mean luminosi- 
ties of the two samples are quite similar: for the Z03 
sample, (Lr) = 1.7 x lO^^/i^^ig, while the mean for our 
sample is (Lr) = 1.46 x 10^°ft L©, both comparable to 
= 1.54 X 1O^°/i-2L0 (Blanton et al. 2001). Assum- 
ing that ro for ^gg scales with luminosity as in Norberg 
et al. (2001), we expect the auto-correlation length for 
our sample of galaxies to be only slightly lower than that 
of Z03, ro = 5.60 rather than 5.77/1"^ Mpc. 

Figure 12 shows ^gm (points with error bars) and the 
^gg from Z03 (solid line) . The S,gg was inferred by invert- 
ing from ujPg using the same techniques outlined in §2.1 
(I. Zehavi, 2003, private communication). The ampli- 
tude of ^gg has been lowered by about 5% to account for 
the slight difference in mean luminosity for the two sam- 
ples as discussed above, ^gg and ^gm for these samples 
are consistent within the errors over the common range 
of radii. The inverse galaxy-mass bias r/b = ^gm/^gg is 
shown in the bottom panel of figure 12. We use the in- 
verse because ^gm is much noisier than ^gg. To calculate 
b/r, we have interpolated the points and errors from Z03. 
The bias b/r is consistent with unity and is scale inde- 
pendent within the errors. Fitting the bias to a constant 
over the radial range 0.2 — 6.7h~^ Mpc, we find 

(6/r) = (1.3±0.2)(^^). (32) 

7. GALAXY-MASS CORRELATIONS IN A VOLUME 
LIMITED SAMPLE OF LUMINOUS GALAXIES 

In the previous section we used a flux limited sample 
to measure ^gm- This selection has the advantage that it 
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Fig. 11. — Correlation matrix for ^gm in figure 10, propagated 
from the covariance matrix of AS. 
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Fig. 12. — Galaxy-mass correlation function (filled circles) and 
the galaxy auto-corrclation function (solid line) for similarly se- 
lected samijlcs of SDSS galaxies. The top panel shows the two cor- 
relation functions and the bottom panel shows the inverse galaxy- 
mass bias rjh = l^gmji^gg- Wc plot the inverse because is 
much noisier than S^gg. The bias is consistent with scale indepen- 
dent, with a mean of ^jr) = (1.3 ± 0.2)(Qm/0.27), shown as the 
hashed region. 



uses the largest possible sample of lens galaxies and thus 
provides the best possible S/N measurement of ^gm- The 
interpretation is somewhat complicated, however, by the 
selection function. In order to provide a measurement 
that is simpler to interpret, we also defined a volume 
and magnitude limited sample, using the limits for the 
high- luminosity sample of Zehavi et al. (2002): 0.1 < 
z < 0.174, -23.0 < Mr - 51og(/i) < -21.5, (L^) = 
3.96 X 1O^°L0. This sample contains 10,277 lenses. 

The inferred ^g„i for this sample is shown in the top 
panel of figure 13, along with S^gg inferred from the same 
volume and magnitude limited sample from Zehavi et al. 
(2002). Power law fits to are shown in table 1. The 
slope is significantly steeper than for the overall sample. 
^gm shows deviations from a pure power-law, however, 
especially at small separations (r < 0.2h^^ Mpc), and 
this is reflected in the poor reduced for the power 



law fit. We will study the shape of ^gm as a function of 
luminosity in more detail in §8, where we split the overall 
sample into multiple bins in luminosity. 

The bottom panel of figure 13 shows the inverse bias 
r/b over the common range of scales: 0.4 < r < 6.7h~^ 
Mpc. The b/r for this volume limited sample is consistent 
with scale-independent. The mean bias over these scales 
is 

{b/r)vum = {2.0 ±0.7) (^^^ (33) 
Because this is a higher luminosity sample than the over- 



all flux limited sample, one expects a larger bias (Norberg 
et al. 2001), but we do not yet detect any significant in- 
crease in the bias, as {b/r)vum is consistent with {b/r) 
at the 1-C7 level. 

8. DEPENDENCE OF GALAXY-MASS CORRELATIONS ON 

LUMINOSITY 

To study galaxy-mass clustering as a function of galaxy 
luminosity, we split the lens sample into three bins of lu- 
minosity separately in each of the five SDSS bandpasses, 
as shown in table 2. The binning is chosen primarily to 
yield comparable lensing S/N in each bin. To achieve 
this we place 85% of the galaxies in the lowest luminos- 
ity bin, 10% in the middle, and 5% in the highest bin. 
The low-luminosity bin must be so large because intrin- 
sically faint galaxies are weakly clustered and are only 
detected in the SDSS at low redshift (for which the lens- 
ing efficiency is low). Since it contains the majority of 
lens galaxies, the faintest subsample has a mean luminos- 
ity comparable to the overall sample (~ L*), while the 
most luminous subsample is ~ 5 — 6 times brighter. Note 
that the number in each bin varies as a function of band- 
pass because galaxies arc not detected in every bandpass. 
Also note that, while the different bins in luminosity are 
independent, the subsamples in the different bandpasses 
are simply re-samplings of the same lenses. 

The mean AE for each luminosity subsample is shown 
in Figure 14. The profiles have been re-binned from 18 
to 9 radial bins for clarity of presentation. Since the 
lowest luminosity subsample contains the large majority 
of lenses, the signal is comparable to that of the overall 
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Fig. 13. — Same as figure 12, but for a volume and magnitude 
limited sample (0.1 < z < 0.174), (-23.0 < Mr - 51og(/i) < 
—21.5). The solid line is £,gg for the same selection criteria from 
Zehavi et al. (2002). We detect no scale dependence in the bias 
b/r. The mean bias over the common range of scales (0.4 — %.7h~^ 
Mpc) is {b/r)viim = (2.0±0.7)(Q„/0.27). 
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sample in both amplitude and shape. The higher lumi- 
nosity subsamples show a significantly steeper slope on 
intermediate to large scales, and a flattening at smaller 
scales. This trend is seen in each of the bandpasses. 

The AE for each luminosity bin was inverted to 
in the same way as for the overall sample. The result- 
ing S,gm for each subsample is shown in figure 15. Note 
that the shape and relative amplitude of these data are 
independent of the assumed flm- 

^gm for the highest luminosity bins deviate significantly 
from power laws, although this trend is less clear than 
for AS. S^gm at large radius (r > Ih^^ Mpc) is consistent 
across the different luminosity bins, while at intermedi- 
ate radii (0.1 < r < lh~^ Mpc) the amplitude of the 
correlations increases with luminosity. On the smallest 
scales (r < 0.1h~^ Mpc), the clustering increases weakly 
with luminosity. These trends are in qualitative agree- 
ment with the predictions from semi-analytic and n-body 
results from Guzik & Seljak (2001). Note, the low data 
points in for the middle (and highest) w-band lumi- 
nosity bin is due to a single negative point in the un- 
binned AS. 

We fit the £^gm to power law models £,gm = {r/ro)~'^ 
for each of the luminosity subsamples. The results for 
these fits, summarized in table 2, indicate that ^gm for 
the low luminosity samples is a good fit to a power law 
and consistent with the overall sample. The two high- 
est luminosity bins have a steeper logarithmic slope than 
the low luminosity sample. The highest luminosity bin, 
however, is a poor fit to a power law in each of the band- 
passes. The fits over all radius give a slope of 7 ~ 2. 
The slope for r > 0.1 Mpc is 7 ^ 2.5, much steeper 
than the slope including the points at small radii. This 
is indicative of the rollover seen at small separations. 

A priori, one might hypothesize that the change in 
slope from to brighter galaxies is partly due to the 
fact that the higher luminosity subsamples are also red- 
der: for the r-band samples, g — r = 0.73(0.77) for the 
middle (brightest) sample, while 5 — r = 0.63 for the full 
sample (see table 2). However, the same increase of slope 
is also seen for the u subsamples, for which the trend of 
g — r color with luminosity is much less pronounced. 

9. DEPENDENCE OF GALAXY-MASS CORRELATIONS ON 
GALAXY SPECTRAL TYPE AND COLOR 

In addition to luminosity, we can also study how the 

lensing signal varies with galaxy spectral type and color. 
The SPECTRDld pipeline classifies galaxies as early to 
late spectral types using the method developed in Con- 
nolly et al. (1995) and Connolly & Szalay (1999). A 
Karhunen-Loeve decomposition (KL) is performed on a 
large set 100, 000) of galaxy spectra, ranging from 
quiescent E-Sc types to starburst irregulars. The spec- 
trum of each galaxy is expanded in terms of the KL 
eigenspectra, and the first five coefficients of the ex- 
pansion {ECOEFFl, . . ,EC0EFF5} are measured. Since 
most of the information is contained in the first two 
coefficients, a simple one-parameter family, constructed 
from the angle in the plane of the first two eigen- 
vectors, provides a sufficient classification: ECLASS = 
tan~^(-EC0EFF2/EC0EFFl). Early types have more neg- 
ative ECLASS. 

The distribution of ECLASS for the lens galaxy sample 
is shown in the bottom panel of figure 16. The rest frame 



g — r distribution is also shown in the top panel. The 
g — r color and ECLASS correlate well. We make cuts at 
ECLASS=-0.06 and g — r = 0.7 to split the sample into 
early/late types and red/blue. These cuts each divide 
the full sample roughly in half. 

The mean AS for these different classes is shown in the 
top panel of figure 17. The red galaxies show a power 
law AS for R > 100 kpc and a fiattening at smaller 
separations, similar to the two highest luminosity bins of 
the full sample. The blue and late type measurements, 
although noisy, are consistent with a power law AS over 
all separations. 

The bottom panel of figure 17 shows the inversion of 
AS to ^grn for each class. The non-power law behavior 
for the early types is still visible for ^gm but it is less pro- 
nounced. The strong low feature in the blue/late samples 
is related to the feature at smaller radius in AS; recall 
that the inversion to ^g„i involves the derivative of AS. 

We fit a power law to each sample, using the full 
covariance matrix for the red/early samples, but only 
the diagonal elements for the blue/late samples because 
their covariance matrix is quite noisy. The results are 
shown in Table 1. The red/early type galaxies have a 
higher amplitude than the blue/late types, as indicated 
by the best-fit correlation length vq. The red/early sub- 
samples are also nearly twice as luminous as the blue/late 
siibsamples. We detect no appreciable difference in the 
slope of the blue/late correlation function relative to that 
of the red/early types, in contrast to the strong shift in 
slope of the galaxy auto-correlation function as a func- 
tion of type and color (Zehavi et al. 2002). On the other 
hand, the errors in the blue/late sample slopes are siz- 
able, reflecting the lower S/N. 

9.1. Dependence of S^gm on Luminosity for Red Galaxies 

The lensing S/N for the red galaxy sample is quite 
high, allowing us to further subdivide this sample into 
three bins of luminosity in each of the 5 SDSS band- 
passes, as we did in §8. The split is again 85%, 10%, and 
5% of the galaxies in each of the three luminosity bins. 
The mean AS for these luminosity bins is shown in Fig- 
ure 18, and the inversion to $_gm is shown in figure 19. 
The trends of AS and ^gm with red galaxy luminosity 
are similar to those for the full sample. The fits to power 
law correlation functions are listed in Table 3. Note that 
the correlation length ro and mean luminosity for each 
red luminosity subsample is larger than that of the cor- 
responding full luminosity subsample, as expected from 
§9. 

We now have two pieces of evidence suggesting that the 
increase in clustering strength and change in shape we 

observe is more closely associated with luminosity than 
with color: (i) for the subsamples split only by luminos- 
ity (not color), galaxies split by ?i-band luminosity show 
a significant increase of clustering strength with luminos- 
ity, and a change in shape, even though the mean g — r 
color is nearly independent of luminosity for these sub- 
samples; (ii) within the red galaxy sample itself, we find 
the same trends that we see for the full sample. 

9.2. Dependence of £,gm on Velocity Dispersion 

For a large fraction of the early type galaxies, se- 
lected by ECLASS and by the requirement that the sur- 
face brightness approximately fits a de Vaucouleurs pro- 
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Fig. 14. — Mean AS in three luminosity subsamples for each of the 5 SDSS bandpasses. In each panel, circles connected by solid 
lines(black), crosses connected by dotted lines(blue), and triangles connected by dot-dashed lines(red) represent measurements for the 
lowest, middle, and highest luminosity subsamples. The data have been re-binned to 9 bins from 18 for clarity. 
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Fig. 15. — Same as figure 14, but now plotting ^gm obtained by inverting AS. The highest luminosity bin shows significant deviation 
from a power law. 
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Fig. 16. — Distribution of g — r (top panel) and ECLASS (bottom 
panel) for lens galaxies. Early types have more negative ECLASS, 
which has been plotted in reverse for comparison with g — r. We 
divide the samples at g — r = 0.7 and ECLASS=-0.06 as shown by 
the gray vertical lines. 



file, the SDSS spectroscopic pipeline measures the veloc- 
ity dispersion ay The distribution of velocity disper- 
sions is shown in figure 20. We use only galaxies with 
50 < fit, (km/sec) < 400 for this study and split the sam- 
ple roughly in half at cr„ = 182 km/sec. The 1ow-ct^ 
sample has a mean of ay — 143 km/s and the high-cr^ 
sample has a mean of = 215 km/s. 

In Figure 21 we show AE and ^gm for these two sam- 
ples of early-type galaxies split by velocity dispersion. 
The half of the sample with higher shows a steeper 
slope, and a larger amplitude at r < 500h~^ kpc, than 
the low-cr„ galaxies. Interestingly, the amplitude for the 
low ay galaxies is actually somewhat higher that that for 
high-fft, galaxies for r > IMpc. For neither sample is 
Cgm well fit by a power law. The shape of the high-Uu 
£,gm is similar to that for the high luminosity red subsam- 
ple. Given the Faber- Jackson relation between velocity 
dispersion and luminosity, this trend agrees qualitatively 
with the results above for the luminosity-scaling of £,gm 
for red galaxies. 

10. DISCUSSION 

We have presented measurements of the mean lensing 
signal AE = E(< R) — E(i?) for a flux limited spec- 
troscopic sample of SDSS galaxies with (L) ^ L* and 
(z) ~ 0.1 over scales 25h~^ kpc to 10h~^ Mpc. AT, is 
consistent with a power-law for all separations, with in- 
dex 0.76 ± 0.05. 

We performed a direct inversion to the three- 
dimensional galaxy-mass cross-correlation function 
over scales 0.025— 6. 7/i~^ Mpc^. The resulting is con- 
sistent with a power-law galaxy-mass correlation func- 
tion over this range, ^gm x (17^/0.27) = (r/ro)~^. The 
measured slope, 7 = 1.79 ± 0.06 is consistent with the 
slope of the auto-correlation function ^gg for a similarly 
selected set of galaxies (Zehavi et al. 2003). We find that 

^ Modelers should predict AS if possible. It is higher S/N and 
contains more information than the inversion. 
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Fig. 17. — Mean AE and ggm for different classes of galaxies. 
The top panel shows AS for galaxies classified as early or red in 
thick solid and dashed lines respectively (red). Late type and blue 
galaxies are shown in thin solid and dashed lines (blue). The trends 
for the ECLASS and g — r cuts are very similar, as expected since 
color and spectral type are highly correlated. The bottom panel 
shows the inversion to ^gm- 



the measured scale length, tq = (5.4±0.7)(0.27/f}™)i/T, 
is also consistent with that of £,gg for the fiducial value of 
rim = 0.27, although the range of radii probed by the two 
data sets is quite different: the galaxy clustering results 
are measured for r > 0.2h~^ Mpc. 

We compared £,gm with £^gg from Zehavi et al. (2003) 
over the common range of radii: 0.2 — 6.7h~^ Mpc. The 
ratio of these functions is S,gg/^gm — b/r, where b is the 
bias and r is the correlation coefhcient at a fixed sepa- 
ration. We find that b/r is consistent with being scale 
independent over this range of scales, with a mean value 
of (b/r) = (1.3 ± 0.2)(f7,„/0.27). Semi-analytic galaxy 
formation models suggest that b/r is a direct measure 
of the standard bias &, i.e., that r approaches unity for 
separations larger than 1 Mpc (Guzik & Seljak 2001), al- 
though there are hints from other lensing data that this 
may not be the case (Hoekstra et al. 2002a). 

We repeated this analysis for a volume and magnitude 
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Fig. 19. — Same as figure 15 but only for red galaxies, g — r > 0.7. 
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limited sample of 10,277 luminous galaxies, with selec- 
tion criteria designed to match that of the most lumi- 
nous sample of Zehavi et al. (2002): 0.1 < z < 0.174, 
-23.0 < Mr - 51og(ft,) < -21.5. We detect no scale de- 
pendence in b/r, although ^g„, for this sample is nois- 
ier than the overall sample. We find a mean h/r of 
(2.0 ± 0.7) (rim/0.27) over the common range of scales 
(0.4-6.7/1-1 Mpc). 

These results are consistent with and extend previous 
galaxy lensing measurements. Fischer et al. (2000) found 
a power law AS for angular scales 9 < 600", or roughly 
R < Ih^^ Mpc. Wilson et al. (2001) and Hoekstra et al. 
(2003a) (which incorporates some lens redshifts) found 
similar results over a smaller range of scales, although 
the power law index was less well constrained. The 
high precision measurements of Hoekstra et al. (2003b) 
show a power law on ~2 Mpc scales (no power law fit 
is quoted). McKay et al. (2002), which used a smaller 
sample of SDSS lens galaxies with known redshifts, found 
a power law AS with slope consistent with the current 
study for 0.025 < R < lh~^ Mpc. An model-dependent 
SDSS lensing analysis concentrating on small scales (Sel- 
jak et al. 2004, in preparation) found results consistent 
with those presented here. Hoekstra et al. (2002a) used 
a deeper, narrower sample than ours, without lens red- 
shifts, to obtain a high S/N measurement of the shear 
over effective scales R < 5h~^ Mpc. They find that b/r is 
approximately scale-dependent over the range 0.1 — 3h~^ 
Mpc, with mean amplitude b/r — 1.09 ± 0.04 over these 
scales. However, combining their galaxy lensing data 
with deeper cosmic shear data, they measured both b 
and r. The errors in b and r are correlated, since both 
depend on the cosmic shear, but they found that b and 
r each vary with scale and track one another in such a 
way that the ratio b/r is approximately constant. 

The value of b/r inferred above requires only the fol- 
lowing assumptions: that general relativity is valid, the 
cosmological principle holds, and the FRW metric prop- 
erly relates redshift to distance in our universe. We have 
explicitly shown the dependence on the parameter 
(the dependence on other parameters through the angu- 
lar diameter distance is weak for these redshifts). Further 
interpretation of these measurements requires further as- 
sumptions and model dependence. For example, if one 
believes that galaxies are nearly unbiased tracers of 
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Fig. 20. — Histogram of measured velocity dispersion (7„ for 
49,024 early-type lens galaxies. 
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Fig. 21. — Mean AS and ^gm for early type galaxies with 
measured velocity dispersion t7„ > 182 km/sec (red-dashed curve) 
and < 182 km/sec (blue-solid curve). The amplitude for high ay 
galaxies is stronger at small scales, as expected from the trend with 
luminosity. 



the mass on these scales, say b/r = bun < 2, then one 
can place a limit flm < 0.6 at 95% confidence. We leave 
it to the reader to further interpret these measurements 
in the context of particular models. 

Splitting our lens galaxies into subsamples has enabled 
us to study £,gm as a function of luminosity, color, and 
spectral type. The amplitude of increases with in- 
creasing luminosity on intermediate (0.1 — Ih^^ Mpc) 
scales, consistent with the lower S/N results of McKay 
et al. (2002). However, we find that S^gm is nearly inde- 
pendent of luminosity on scales larger than lh~^ Mpc 
(see figure 15). This is indicative of the increase in slope 
with luminosity on intermediate scales. The highest lu- 
minosity bin shows significant deviation from a power 
law, with a fiattening at radii less than 0.1h~^ Mpc. 
These trends of the slope and amplitude of ^gm with lu- 
minosity are in good qualitative agreement with the pre- 
dictions of semi-analytic galaxy formation models (Guzik 
& Seljak 2001). The semi-analytic models predict that 
the amplitude of increases with luminosity on small 
scales, but is only weakly dependent on luminosity for 
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scales r > Ih^^ Mpc, in agreement with our resuhs. 

The increase in clustering amplitude with luminosity 
on intermediate scales for is similar to that seen for 
the galaxy-autocorrelation function (Norberg et al. 2001; 
Zehavi et al. 2002), but the deviation from a power law 
on small scales at higher luminosity has not been seen in 
previous autocorrelation measurements. Also, the auto- 
correlation function shows a nearly equal increase in am- 
plitude with luminosity on both large and intermediate 
scales (Norberg et al. 2001; Zehavi et al. 2002), but ^g^, 
does not. This may imply a scale-dependent bias for high 
luminosities. 

We have also detected a trend in with galaxy color 
and spectral type (figure 17). Red/early type galaxies 
show stronger correlation amplitude than blue/late types 
but the slopes of £_g„i for the two subsamples are consis- 
tent within the errors. These results are consistent with 
the previous results of McKay et al. (2002). This behav- 
ior is also in qualitative agreement with the semi-analytic 
galaxy formation predictions (Guzik & Seljak 2001). 

Finally, we have sufficient S/N to study the scaling of 
Cgm with luminosity and velocity dispersion for red/early 
type galaxies. Again, we find a similar increase in loga- 
rithmic slope and amplitude with luminosity on interme- 
diate scales, and a weak luminosity dependence on large 
scales. We find the slope of ^gm increases with cr„, but 
the high sample shows significant deviations from a 
power law. 

A separate analysis will compare these results in de- 
tail with predictions from N-body simulations and the 
halo model of structure formation. A future study, using 
a larger lens sample, will concentrate on galaxy-galaxy 
lensing at small separations, where the signal is primarily 
sensitive to galaxy dark matter halos. 
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11. APPENDIX: PHOTOZ SYSTEMATIC ERRORS AND 
INTRINSIC ALIGNMENTS 

We discussed two sources of systematic errors in §4.1: 
residuals from the PSF correction and clustering of faint 
sources around the lens galaxies. These biases have been 



estimated from the data and corrected for in our final 
results. 

Another possible source of error may come from biases 
in the photometric redshifts (§3.1.5). Although the pho- 
tometric redshift distributions inferred for SDSS galaxies 
are consistent with published redshift surveys for r < 21 
(Csabai et al. 2003), there are well-known degeneracies 
in the technique (between type and redshift) that tend to 
push galaxies to particular values of photoz. This pile-up 
can be seen in figure 2: there are significant peaks in the 
photoz distribution at several redshifts. Although such 
peaks are often seen in spectroscopic redshift surveys, 
where they are associated with large-scale structures, the 
volume probed by the SDSS photoz sample is so large, 
containing ~ 4 million galaxies covering ~ 1000 square 
degrees out to 2; ~ 0.6, that the effects of such structures 
should be very small. 

Because galaxies in parameter regions with strong de- 
generacies have relatively large photoz errors, we expect 
this bias to be suppressed in the final analysis. We have 
repeated our analysis using only the source galaxy red- 
shift prior inferred from summing the individual Gaus- 
sians (as shown in figure 2), i.e., not using the photoz 
measurements for individual source galaxies, and we re- 
cover the same results within the noise. We have also 
used the method outlined in McKay et al. (2002), which 
uses the r magnitudes to infer the redshift distribution, 
and recover consistent results. Thus the bias from pho- 
toz degeneracies appears to be small compared to our 
statistical uncertainty. 

Another issue that can possibly complicate interpre- 
tation of our results is intrinsic correlations in the el- 
lipticities of the source galaxies. If the orientations of 
physically associated pairs of galax:ies are correlated, this 
could possibly mimic a lensing signal. 

Tidal torques can cause the angular momentum vec- 
tors of physically associated galaxies to be aligned (Crit- 
tenden et al. 2001). The magnitude of these ellipticity 
correlations is expected to be strong and relatively con- 
stant on scales less than ~ 1 Mpc, and decrease as ^ls, 
the correlation function of sources with lenses, on large 
scales (Crittenden et al. 2001). This effect has been seen 
in N-body simulations (Croft & Metzler 2000; Heavens 
et al. 2000) and there are possible detections in the TuUy 
catalog (Pen et al. 2000) and SuperCOSMOS Sky Sur- 
vey data (Brown et al. 2002), although these data are 
not ideal for such a study. The recent work of Heymans 
et al. (2003) suggests this effect may be near the low end 
of theoretical predictions. 

It is not clear how intrinsic alignments affect tangen- 
tial shear measurements. The tangential shear is the 
correlation between the position of a lens and the pro- 
jection of source ellipticities onto the tangential reference 
frame, defined by a circle centered at the lens. The in- 
trisic ellipticity correlations described above are correla- 
tions between pairs of ellipticities. There is little theoret- 
ical guidance about what tangential alignment to expect. 
For limits on this effect we turn to the recent study by 
Bernstein & Norberg (2002) which uses data from the 
2dF redshift survey to study the tangential alignment of 
faint (2.2 magnitudes fainter than the primary) physi- 
cally associated galaxies. They place an upper limit of 
(e!p*) < 0.02(95%) within 500 /i'^ kpc at (z) = 0.1. A 
recent analysis of this tangential alignment within the 
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SDSS spectroscopic sample (Koester et al. 2003) places 
still more stringent limits on this tangential alignment 
(e!p*) < 0.007(95%) on 500/1"! y.^^ 

Most of the sources used in our study are far behind the 
lenses, in which case no correlation is expected between 
lens position and the tangential source galaxy ellipticity. 
A small fraction of sources, however, arc actually faint 
companions of the foreground lenses. The fraction of 
physically associated sources is generally small (Fischer 
et al. 2000), and the use of photometric redshifts in our 
study down- weights these sources in the final analysis. 
The relative contribution of these objects in the Icnsing 
measurement is shown in figure 9. The mean contribu- 
tion over 500 kpc is 0.018. 



Assuming the result from Bernstein & Norberg (2002) 
also applies to our faint sources, wc place a limit on the 
contamination from intrinsic alignments at (e!p*) < 0.02* 
0.018 = 3.6 X 10"^ on 500/1"^ kpc scales. The mean 
distortion within this radius is 1.5 ± 0.1 x 10~^. Thus, 
on 500 /i~! kpc scales, the intrinsic alignments should 
contribute less than 20% of the signal at 95% confidence. 
The new SDSS results suggest the eflFect is smaller still (< 
8.5%). The intrinsic alignment correlations decrease with 
scale roughly as and the relative contribution from 
the sources also decreases quickly with scale (figure 9). 
Thus, the contribution from intrinsic alignments should 
be small for radii r > 500h~^ kpc. 
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Table 1. Model Fits for Different Lens Samples 



Sample 


Selection Criteria 


Mean Abs. Mag. 


Mean g — r 


^Lenses 


ro 


7 




All 




-20.767 (1.455 ± 0.004) 


0.629 


127001 


5.4 ± 0.7 


1.79 ± 0.06 


17.1/15 


Red 


g-r > 0.7 


-21.061 (1.908 ± 0.006) 


0.753 


60099 


6.9 ± 0.8 


1.81 ± 0.05 


12.3/15 


Blue 


g-r < 0.7 


-20.477 (1.114 ± 0.004) 


0.536 


65134 


4.0 ± 1.0 


1.76 ± 0.15 


5.97/6 


Early 


ECLASS < -0.06 


-21.036 (1.864 ± 0.006) 


0.737 


62340 


7.3 ± 0.8 


1.77 ± 0.05 


16.9/15 


Late 


ECLASS > -0.06 


-20.474 (1.111 ± 0.004) 


0.539 


64378 


3.3 ± 1.0 


1.89 ± 0.17 


7.18/6 


Vlim 


-23.0 < Mr < -21.5 


-21.854 (3.961 ± 0.012) 


0.718 


10277 


5.1 ± 0.8 


2.01 ± 0.07 


24.6/15 




O.K z < 0.174 













Note. — Absolute magnitudes are r-band Pctrosian M — 5 logjQ h. Values in parentheses are luminosity in units of 10^'^ h ^Lq . 
The means are calculated using the same weights as the lensing measurement. The value of Mt,{Lt) for the r-band is -20.83 
(1.54). The ro and 7 are best fit parameters for ^gm = (r/ro)~^; ro is measured in Mpc. A value of Qm. = 0.27 was assumed. 
For the "late" and "blue" samples, the data were rebinned from 17 to 8 radial bins. 



Table 2. Luminosity Bins for All Galaxies 



Bandpass Abs. Mag. Range Mean Abs. Mag. Mean g — r N^enses ro 7 X^/'^ 



u 


-19.6 


< 




< 


-15.0 


-18 


515 


(0 


908 


± 


0.002) 


0.628 


106585 


5.5 


± 


0.8 


1 


75 


± 


0.06 


13.5/15 




-20.0 


< 


M„ 


< 


-19.6 


-19 


763 


(2 


864 


± 


0.003) 


0.639 


12539 


4.4 


± 


0.8 


2 


04 


± 


0.09 


19.2/15 




-22.0 


< 


Mu 


< 


-20.0 


-20 


318 


(4 


779 


± 


0.020) 


0.654 


6270 


6.9 


± 


1.0 


1 


97 


± 


0.07 


26.3/15 


9 


-21.0 


< 


Mg 


< 


-16.5 


-19 


891 


(0 


956 


± 


0.002) 


0.623 


106646 


5.4 


± 


0.9 


1 


75 


± 


0.07 


15.1/15 




-21.4 


< 


Mg 


< 


-21.0 


-21 


188 


(3 


156 


± 


0.004) 


0.699 


12546 


5.0 


± 


0.8 


2 


02 


± 


0.08 


12.7/15 




-23.5 


< 


Mg 


< 


-21.4 


-21 


718 


(5 


145 


± 


0.025) 


0.734 


6275 


7.0 


± 


0.8 


2 


03 


± 


0.06 


24.6/15 


r 


-21.7 


< 


Mr 


< 


-17.0 


-20 


544 


(1 


184 


± 


0.003) 


0.621 


106643 


5.4 


± 


0.8 


1 


74 


± 


0.07 


15.9/15 




-22.2 


< 


Mr 


< 


-21.7 


-21 


902 


(4 


137 


± 


0.005) 


0.726 


12543 


5.2 


± 


0.8 


2 


03 


± 


0.07 


13.5/15 




-24.0 


< 


Mr 


< 


-22.2 


-22 


463 


(6 


937 


± 


0.026) 


0.767 


6274 


7.5 


± 


0.8 


2 


01 


± 


0.05 


30.9/15 


i 


-22.0 


< 


M, 


< 


-17.0 


-20 


870 


(1 


446 


± 


0.003) 


0.620 


106625 


5.3 


± 


0.9 


1 


74 


± 


0.07 


17.6/15 




-22.5 


< 


M, 


< 


-22.0 


-22 


232 


(5 


067 




0.006) 


0.734 


12544 


5.1 


± 


0.7 


2 


04 


± 


0.07 


19.3/15 




-24.0 


< 


M, 


< 


-22.5 


-22 


776 


(8 


368 


± 


0.030) 


0.772 


6271 


6.9 


± 


0.8 


2 


06 


± 


0.06 


24.7/15 


z 


-22.2 


< 


M^ 


< 


-17.0 


-21 


069 


(1 


721 


± 


0.004) 


0.619 


105975 


5.3 


± 


0.9 


1 


73 


± 


0.07 


16.8/15 




-22.6 


< 




< 


-22.2 


-22 


350 


(5 


599 


± 


0.007) 


0.725 


12466 


5.5 


± 


0.8 


1 


97 


± 


0.07 


20.3/15 




-24.0 


< 




< 


-22.6 


-22 


881 


(9 


132 


± 


0.032) 


0.753 


6236 


6.6 


± 


0.8 


2 


05 


± 


0.06 


22.1/15 



Note. — Galaxies were split into three bins of absolute magnitude in each of the five SDSS bandpasses, as listed in column two. 
See table 1 for explanations of the other columns. The value of M.(L.) is -18.34(0.77), -20.04(1.10), -20.83(1.54), -21.26(2.07), 
-21.55(2.68) for -u, gr, r, i, z respectively. 



Table 3. Luminosity Bins for Red Galaxies 



Bandpass 


Abs. Mag. Range 


Mean Abs. 


Mag. 


Mean g — r 


^Lenses 




ro 










7 






u 


-19.6 < Mu < 


-15.0 


-18.521 


(0.913 


± 


0.003) 


0.752 


51074 


6.9 


± 





9 


1 


77 


± 


0.06 


10.0/15 




-20.1 <Mu< 


-19.6 


-19.791 


(2.939 


± 


0.005) 


0.769 


6008 


5.6 


± 





8 


2 


09 


± 


0.07 


21.6/15 




-22.0 <Mu< 


-20.1 


-20.373 


(5.026 


± 


0.027) 


0.789 


3006 


9.6 


± 


1 


1 


1 


96 


± 


0.06 


35.2/15 


9 


-21.2 < Mg < 


-16.5 


-20.112 


(1.171 


± 


0.003) 


0.752 


51084 


7.0 


± 


1 





1 


75 


± 


0.06 


9.46/15 




-21.6 <Mg< 


-21.2 


-21.364 


(3.712 


± 


0.006) 


0.773 


6008 


6.4 


± 





7 


2 


04 




0.06 


20.8/15 




-23.5 <Mg < 


-21.6 


-21.841 


(5.762 


± 


0.025) 


0.788 


3006 


8.9 


± 





9 


2 


03 


± 


0.05 


40.2/15 


r 


-21.9 < Mr < 


-17.0 


-20.867 


(1.596 


± 


0.005) 


0.751 


51080 


7.0 


± 


1 





1 


75 


± 


0.06 


9.56/15 




-22.4 < Mr < 


-21.9 


-22.140 


(5.152 


± 


0.009) 


0.775 


6008 


6.6 


± 





8 


2 


03 


± 


0.06 


19.0/15 




-24.0 < Mr < 


-22.4 


-22.635 


(8.127 


± 


0.037) 


0.793 


3004 


9.0 


± 





9 


2 


03 


± 


0.05 


34.3/15 


i 


-22.3 < Mj < 


-17.0 


-21.240 


(2.033 


± 


0.006) 


0.751 


51078 


6.9 


± 


1 





1 


76 


± 


0.06 


10.1/15 




-22.7 < < 


-22.3 


-22.478 


(6.354 


± 


0.011) 


0.776 


6007 


7.1 


± 





9 


1 


98 


± 


0.06 


20.9/15 




-24.0 <Mi < 


-22.7 


-22.946 


(9.778 


± 


0.043) 


0.790 


3005 


8.4 


± 





8 


2 


06 


± 


0.05 


35.0/15 


z 


-22.4 <Mz< 


-17.0 


-21.412 


(2.360 


± 


0.007) 


0.751 


50697 


7.0 


± 


1 





1 


75 


± 


0.06 


10.1/15 




-22.8 < M^ < 


-22.4 


-22.563 


(6.810 


± 


0.011) 


0.772 


5962 


5.4 


± 





9 


2 


06 


± 


0.08 


8.99/15 




-24.0 < M^ < 


-22.8 


-23.045 


(10.61 


± 


0.048) 


0.782 


2983 


7.2 


± 





9 


2 


09 


± 


0.06 


26.0/15 



Note. — 



Same as table 2 but for red galaxies, defined as galaxies with g — r > 0.7. 



